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I. INTRODUCTION 

Over the last two decades, solar neutrino physics evolved from a qualitative description of the data into a precision 
science. Following the seminal experiment of Davis [1] , experiments were performed to measure the predominant pp 
neutrinos [21 [3], as well as the higher energy neutrinos: Borexino measured the line flux coming from ^Be neutrinos 
[3]. The continuous flux of the ^B neutrinos were measured at SuperKamiokande (SuperK) |5j and at Sudbury 
Neutrino Observatory (SNO) [B| with a recent confirmation from the KamLAND experiment [7]- The Cerenkov 
detector experiments SuperK and SNO are high statistics experiments. In addition, SNO experiment was successful 
in measuring the total high-energy solar neutrino flux using two different methods to detect the neutral current breakup 
of deuteron [8]. This measurement, providing the total flux independent of neutrino mixing and matter effects, is 
in general agreement with the predictions of the Standard Solar Model [51 [TO]. In parallel with these experimental 
developments there have been increasingly refined analyses of the solar neutrino data alone or combined with the 
data from accelerator, reactor, and atmospheric neutrino measurements (see e.g. Refs. [HHHIIIS])- In the particular 
case of solar neutrinos, two distinct regions of the neutrino parameter space are identified: the so-called LOW region 
corresponding to smaller values of Sm^i 10^^ eV^) and the so-called LMA region^ corresponding to larger values 
of Sm2i (~ 10~^ eV^). The LMA region was originally adopted as the appropriate region of the parameters space 
because of the KamLAND reactor neutrino observations. 

The original analyses of the various phases of the SNO experiment were performed with a rather large energy 
threshold. More recently SNO was able to lower its threshold to an effective electron kinetic energy of 3.5 MeV 
[TT] . With such high precision reached in detecting solar neutrino flux and associated observables, it is worthwhile 
to visit various matter effects and their ramifications. Furthermore, our measurements of the solar neutrino flux are 
not complete: neutrinos from the CN cycle [T5] are not yet identifled. (These neutrinos could probe the interaction 
between the early Sun and its protoplanetary disk [TO]). The sun can also be used as a laboratory for searching 
for physics both beyond the standard solar model and the Standard Model of the particle physics. Examples of the 
former includes searching for effects of density fluctuations [201 HI] ■ Examples of the latter include the effects of large 
neutrino magnetic moments [22]; neutrino interactions beyond those in he Standard Model such as flavor-changing 
neutral currents [23], mass- varying neutrinos [24', and new long-range forces [26j . A search for the subdominant 
physics operational in the Sun requires a thorough understanding of the dominant matter and neutrino mixing effects. 
For example as we demonstrate below, the high-statistics solar neutrino experiments SNO and SuperK cannot really 
distinguish between LOW and LMA solutions. Breaking this degeneracy requires using the KamLAND data. To 
completely rule out the unlikely possibility of CPT violation (i.e the possibility of 6m^ being different for solar 
neutrinos and the reactor antineutrinos), one has to examine another matter-effect, namely the day- night asymmetry. 
Indeed Borexino experiment recently measured this asymmetry for the ^Be neutrinos, verifying the LMA solution 
(and supporting the CPT invariance) using only solar neutrinos [37] . 
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^ LMA stands for "large mixing angle". This is somewhat of a misnomer as the LOW region is also operative for large mixing angles. 
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FIG. 1: Solar electron neutrino survival probability as a function of 5rn? with cos26„ = 0.36 and = 10 MeV. The average 
solar electron density at the neutrino production region is taken to be A'^^ ~ 100 Na/cxv? ■ The horizontal line depicts the value 
measured by SNO at the low-threshold analysis [17j . 



In the next section we provide a critical look at the SNO low-threshold analysis and Borexino Day- Night Asymmetry 
measurement. We show that the solar neutrino data constrains the mixing angle poorly. In section III, we examine 
the effects of new physics beyond the Standard Model and show that subdominant Standard Model effects can mimic 
new physics. (Mathematical details are included in an Appendix). Finally, a brief discussion of our results conclude 
the paper. 



II. SNO LOW-THRESHOLD ANALYSIS AND BOREXINO DAY-NIGHT ASYMMETRY 

MEASUREMENT 

Low-energy threshold analysis of the phase I and phase II data sets of the Sudbury Neutrino Observatory was also 
reported as the electron- neutrino survival probability expanded around the neutrino energy = 10 MeV [17] . In 
this section we first show that, with most mixing angles, the solar electron-neutrino survival probability is doubly 
degenerate in Sm^i- Assuming two-flavor mixing, the electron neutrino survival probability detected on the Earth 
during the day and averaged over the Earth-Sun distance is [55] 

P{iye^iye) = ^ + ^COs26'„(cOS26',;)so„^ce (1 - 2|*2,(S)H 

- ^cos26i^,(sin26i,),o„^ce (*i,(5)^2,(s) +**i,(s)*2,(s)) ' (II-l) 

where 9^ and 9i are the vacuum and initial matter mixing angles respectively, (• • • ) source denotes averaging over the 
neutrino production region in the Sun, and 5*1, (s) a-nd 5'2,(S) are the matter eigenstates calculated on the surface 
of the Sun. The quantity |^2.(s)P is sometimes referred to as the hopping probability, Phop- In the regions of the 
neutrino parameter space consistent with the solar neutrino observation one has either cos29i ^ — 1 (the so-called 
LOW region) or 5'2,(S) ~ (the so-called LMA , or large mixing angle, region). In both cases the last term in Eq. 



(II.ll can be dropped and one obtains [29] 



P^iiy^ -^'''^) = \ + \ cos 29y {cos 20,) source (1 - 2Phop) ■ (11.2) 

For a given value of the mixing angle, the energy (or more correctly 5m? jE) dependence of the survival probability 
comes from the product cos26'i(l — 2Phop)- For smaller values of 5w? (^ 10^'' eV^), even though the matter angle 
is at its maximal value, Phop is finite (but small) for a range of neutrino energies. For larger values of Sir? (~ 10~^ 
eV^), Phop vanishes, but the matter angle is no longer maximal, hence slowly changes with neutrino energy. Hence 
for a fixed value of the vacuum mixing angle, there are two values of Sm? to obtain a given solar electron neutrino 
survival probability as depicted in Figure [l] 
The cosine of the matter mixing angle is 

oa {2BE-COS29,) 
cos29m = / (II-3) 

{2BE-cos29y y+sm^ 29y 
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FIG. 2: A comparison of the hopping probabilities calculated for an exponen tial e lectron density. The exact result given by 
Eq. (II. 6l (thick line) is plotted versus the approximate result given by Eq. (II. 7l (thin line). The same neutrino parameters 
are used as in Fig. [l] 



where, for later convenience, we defined 



B = V2GFNe/Sm^. 



The Taylor expansion of Eq. (II. 3) around Eq — 10 MeV gives 

sin^26lo, , , sin'' 26*0 , ,9 

cos 29m = cos 29„ - 2B . ^„ {E - Eq) - GB^^ — ^ cos 26*0 (-B - £^0) 



sin 26„ 



sin-" 29,, 



(11.4) 



(11.5) 



In the equations throughout this section the subscript denotes a quantity calculated at the neutrino energy of 10 
MeV. An approximate expression for the hopping probability can be estimated by exploiting the fact that the solar 
density decreases more or less exponentially, Ne — A^o exp(— r/rp). For such a density profile the hopping probability 
is given by [30] 



-7r(S(l-cos2e„) 



hop 



1 



(11.6) 



where S — r^Snn? /2E. In the Sun, the arguments of the exponentials of Eq. (II.6) are large and sometimes the 
approximate expression 31j 



hop ' 



-7r(5(l-cos 26^) 



{11.7) 



is used. Hopping probabilities form these two expressions are plotted in Fig. [2j One sees that, for the pertinent 
values of 5m?', the difference between the hopping probabilities predicted by these two expressions are very small. For 
the arguments below we will use Eq. (II.7). (Since the solar electron density is not exactly an exponential, it is also 

II. 7| ) to near exponential densities [35] in the analysis presented below). 
, we get 



possible to utilize a generalization of Eq. 
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Low-energy threshold analysis of SNO reports the daytime solar electron neutrino survival probability as |17| 

P^{iy, ^ ly,) =co + ci{E- 10) + C2{E - 10)^, (II.9) 

where energies are measured in MeV. Note that the values of ci and C2 given by SNO are consistent with zero. In 
this regard, including a third order term in Eq. (II.9) would have been impractical as it would have had a large 
uncertainty. From Eqs. (II.2), (II. 5 1 and (II. 8 1 we expect 
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FIG. 3: Chi-square values of our numerical fit to the SNO LTE results with i 
to distinguish LOW and LMA regions with the SNO data set 



0. As this figure illustrates, it is not possible 



Assuming the values cos26y = 0.36, 6m2i = 
1331 . and using the SNO central value of cq 



7.6 X 10~^ eV^ currently recommended by the Particle Data Group 
: 0.3435 , this gives cos20o = —0.87 which corresponds to an average 
initial density of ^ 100 Na/cui^ in the production region of 10 MeV neutrinos, in excellent agreement with the 
Standard Solar Model [21 [TO]- For the LOW solution we get a hopping probabihty of P^^^ ~ 0.065, yielding 6 = 1.36 
or for ^ 8.2 x 10~^ eV^, we get tq = i?0/lO.6, again in good agreement with an exponential fit to the solar 
electron density profile !34| . For the parameter ci we expect 



2ci 



cos 2d„ 



-2B 



'29o 



sin 2e„ 



for LMA 



-Phop'^ogPhop/Eo for LOW 



(11.11) 



It is widely accepted that studies of the matter effects in solar neutrino oscillations over the last decade established 
Sm2i to be positive [T3] yielding a positive value for B. This makes the sign of ci negative for the LMA solution in 
Eq. (11.11 in contrast with the SNO value. On the other hand the LOW solution has the same sign as the SNO value 
(note that Phop < !)• These considerations suggest that SNO data will have a slight preference for the LOW solution. 
We present numerically calculated for the entire SNO data set in Fig. [s] Clearly qualitative arguments given in 
this section are verified. 

SNO low-energy threshold analysis also provides a fit to the day-night asymmetry. For the LMA and LOW regions 
night-time solar electron neutrino survival probability, averaged over the Earth-Sun distance is given by [351 136] : 



1 



P {Ve Ve) ^ 7,+ - {cos 20^) source (1 " 2Phop) COS 26*^ COS 2(6le - Oy) 



(11.12) 



where O,, is the matter mixing angle inside the Earth. In writing Eq. (11.11 1, matter density in the mantle of the 



Earth (through which neutrinos detected at SNO would travel) is assumed to be constant. The necessary expression 
for the day-night asymmetry is 



A 
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N 



P^ 



pN pD ■ 



(IL13) 



Using the central value of the asymmetry extracted by SNO at E^, = 10 MeV, A = 0.0325, and the initial matter 
angle cos20i = —0.87 in the LMA region or the hopping probability Phop = 0.065 in the LOW region, obtained 
above, we find that the degeneracy of the solution with respect to Sni^i survives both with day and night data. Note 
that the asymmetry quoted by SNO is consistent with zero when the statistical and systematic errors are taken into 
account. We also present numerically calculated day-night asymmetry in Fig. [4j Numerical results agree with the 
simple analytical expression, Eq. (11.13) and the qualitative discussion presented above. 

Recently the Borexino experiment reported a measurement of the day-night asymmetry for the E^ = 0.86 MeV 
''Be line neutrinos j27j . We show the calculated value of this asymmetry in Fig. [5j Clearly the reported value of A = 
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FIG. 4: Calculated day-night asymmetry for the SNO data set. The da shed line is for the exact numerical calculation and the 
solid line is given by the approximate analytical expression, Eq. (11.13 1. 



0.001 ± 0.012 (stat) ± 0.007 (syst) completely rules out the LOW region of the neutrino parameter space using only 
the solar neutrino data without any direct reference to the reactor antineutrinos. This experimentally establishes that 
the neutrino and antineutrino mixings are indeed identical. 

So far we have been discussing only two-flavor mixing, which is only correct for 613 = |37| . Already several years 
ago it was pointed out that the tension between solar neutrino experiments and the KamLAND reactor neutrino 
experiment could be explained by a small, but non-zero value of 613 [38]. SNO low-threshold analysis reports a value 
of 6*13 = 7.2+2;0 de grees, albeit with non-Gaussian errors [T7]. A global analysis by the KamLAND collaboration 
indicates sin^ 6'i3 — O.OOO^qqqy at the 79% C.L. 40,. Best current indication of a non-zero value of 6*13 are given 
by the T2K and MINOS experiments. T2K collaboration announced a 2.5a result: 0.03 < sin^ 26'i3 < 0.28 for the 
normal and 0.04 < sin^ 2^13 < 0.34 for the inverted hierarchies, both for Sep = and maximal 623 PTI. MINOS 
collaboration reports < 2 sin^ 2^13 sin^ ^23 < 0.12 for the normal hierarchy and < 2 sin^ 2^13 sin^ ^23 < 0.19 for 




FIG. 5: Calculated day-night asymmetry for the Borexino data. The dashed line is for the exact numerical calculation and the 
solid line is given by the approximate analytical expression, Eq. (11.13 1. 
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FIG. 6: The 90% confidence level regions of the neutrino parameters allowed by only all the solar neutrino experiments for 
different values of ^13 (in radians). 



the inverted hierarchy, again both for 6cp ^ [15]^. Recent global analyses are consistent with these values [T^ I15j . 
Currently three reactor neutrino experiments (Double Chooz [43l, Daya Bay [44], and RENO [45]) are planning to 
directly measure this angle. 

To discuss the effects of ^13 on solar neutrino data many times the formula [TI] [TBI ES] 

Psxsit'e -> = cos"* 6*13 P2x2(j^e ^ i^e with COS^ 6*13) + siu"^ 6*13 (11.14) 

is used. A discussion of the accuracy of this formula was given in Refs. [46^ and ^7]. We will utilize this formula in 
the discussion given below. 

We present 90% confidence level regions of the neutrino parameters allowed only by all the solar neutrino experiments 
for different values of 6*13 in Fig. [6j An examination of the figure reveals that changing the value of 6*13 has very 
little effect on the dm^i and sin0i2 parameter space. Reasons beyond the insensitivity of the solar neutrino data to 
the value of 6*13 can be explicated by investigating constant P{i'e v^) contours in the 9i2 — 9i3 plane. We present 
constant electron neutrino survival probability contour for a 10 MeV solar neutrino with Smi2 = 7.6 X 10-5 eV2 and 
a constant survival probability for a 0.862 MeV neutrino in the 9i2 — O13 plane in Fig. [7] The survival probability 
for the 10 MeV neutrino is set to the value given by SNO in the low-threshold analysis and the probability for the 
0.862 MeV neutrino is set to the value given by Borexino [15] • In R-ef- [13 it was pointed out that the central values 
from low energy neutrino experiments and high energy neutrino experiments together close in on a small range of 
013- However, the large uncertainties from the low energy experiment mean that even though the two isocontours 
intersect, the range of allowable values for is large. We conclude that the current solar neutrino data alone would 
constrain the mixing angle ^13 rather poorly. 



^ One should note that appearance experiments can only measure the combination sin^ SSia sin^ 6231 hence the uncertainty in the measured 
value of 623 needs to be folded in to quote a separate value of Big. 
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FIG. 7: Constant electron neutrino survival probability contours for a fixed value of Smi2 (taken to be 7.6 x 10~^ eV^) in 
the 012 — Ois plane. The solid line is for a 10 MeV neutrino with survival probability set to the value given by SNO in the 
low-threshold analysis. The dashed line is for a 862 keV neutrino with survival probability set to the value given by the 
Borexino experiment. The wide shaded and hatched areas correspond to one-sigma uncertainties quoted by the Borexino and 
SNO experiments, respectively. 



III. STANDARD MODEL VERSUS NEW PHYSICS IN MATTER EFFECTS 



A convenient way to parameterize the effect of new interactions on neutrinos traveling through media is [531 ISOl [51] 



d 



-cos 2012 sin 26*12 
sin 26*12 cos 26*12 
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(III.l) 



where the terms en and ei2 represent contributions from new physics. Here we wish to emphasize that non-zero values 
of these terms do not necessarily imply new physics. (Clearly en — —1 and ei2 = represent neutrino oscillations in 
vacuum). A derivation of the Eq. (11.14 1 is outlined in the Appendix. In solving the neutrino propagation with ^13 7^ 



one solves an effective two- flavor equation where Nf. is multiplied by cos ^13, which corresponds to en = — sin ^13 
and ei2 — 0. The fact that various oscillation effects may be confused with non-standard interactions was already 
emphasized in Ref. [52 . 

We rewrite Eq. (A.7) for a neutral medium as 



with 



3v^G 



pa 



■jrYe sin^ 6w \ f^w J 



log- 



2 + Ye 

6 



where the electron fraction, Yg, for a neutral medium is given as 



Hence inclusion of the loop corrections would appear in the effective two-flavor evolution as 



eii = - sin^ 6*13 + V^p (sin^ 6*13 cos^ 6*23 - 2 sin^ 9 



23) 



and 



ei2 = sin 6*13 sin 26*23 VVue^ 



(111.2) 



(IH.3) 



(111.4) 



(IH.5) 



(III.6) 



Even though Eq. (III.l) includes both the Standard Model physics and physics beyond the Standard Model, it 
still is a convenient parametrization of the subdominant matter effects. In Figure [8] we show our results. Clearly 
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FIG. 8: The 90 confidence regions for 613 = 0.0 (tlie innermost siiaded area), 0.1 (wider siiaded area and 0.2 (widest sliaded 
area). In making this plot 5m|i is taken to be 7.59 x 10~^ eV^ and tan0i2 = 0.68. 



positive values of ei2 will effectively increase coupling between flavors, whereas negative values would decrease it. 
Since experiments indicate that solar electron neutrino flux is much less than expected, coupling between flavors 
cannot be significantly decreased. Conse quent ly relatively large negative values of ei2 (ei2 < — sm29i2/{2BiE), Bi 
being the value of the parameter of Eq. (II. 4 1 where the neutrino is produced) are not permitted as seen in Fig. |8] 



On the other hand positive values of £12, since they increase mixing, are favored'^ regardless of the value of en. One 
also observes that the allowed parameter space is shifted towards negative values of en. The sign of en can be either 
positive or negative: non-zero ^13 contribution is negative whereas Standard Model loop corrections are positive. To 
understand this shift, we consider the matter mixing angle in the presence of these additional terms 

cos2gi2-2g^(l + en) 

cos2^„j = — III. 7 

y/{2BE{l + en) - cos 26*12)2 + (2B£;ei2 + sin26ii2)2 



and its derivative 




2S^[ei2 cos 2012 + (1 + en) sin 2^12] ^^^^ , 



[2BE{1 + en) - cos 26li2]2 + [2BEei2 + sin26li2]2 

in the presence of these additional terms. For the values of 5m^2 ~ 10~^ eV^, the neutrino evolution is adiabatic as 
can be seen by examining the adiabaticity parameter: 

2BE[ei2 cos 29i2 + (1 + en) sin 2^12] ^^^^ 




{[2BE{1 + en) - cos20i2]2 + [2BEei2 + sin20i2]- 



n3/2- 



For values of en such that the condition 2BE{1 + en) ~ cos20i2 is satisfied, an MSW resonance is possible. In this 
case the adiabaticity parameter takes its maximum value of at the resonance: 

2Si?(l + eii)sin0i2 ^^^^ 




(2B£:ei2 +sin26li2)i/2 



which is small in the Sun for 6m ^2 ^ 10 eV . For larger absolute value of en (both positive and negative) clearly the 



adiabaticity parameter of Eq. (III. 9) is even smaller. (This is also true for large values of 612). Resulting adiabatic 



Of course very large values of ei2 will yield an energy-independent survival probability (0.5) in contradiction with the data, comparing 
low- and high-energy neutrinos. 
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FIG. 9: Percentage change in the electron neutrino survival probability for en = 0.1 (solid line), en — 0.05 (long-dashed line), 
and en = 002 (dotted line). In making this plot (5m|i is taken to be 7.59 x 10~^ eV^, tanSi2 — 0.68, and ei2 = 0. 



neutrino evolution is then controlled by the initial value of the cosine of the matter angle given in Eq. (III.7I. 
Regardless of the value of ei2, for large positive values of en this angle tends to cos29„i ^ —1, whereas for negative 
values of en this cosine may significantly move away from —1. In the absence of the e parameters cos20„j is already 
very close to —1. Consequently more of the parameter space for the negative values of en is available for fitting the 
data, shifting the allowed region of the en towards negative values for a fixed value of ei2. This is in concordance with 
the numerical calculations shown in Fig. [Sj One should point out that values of en and ei2 beyond a few percent are 
likely to be unphysical, but we show relatively large values in Fig. |8]to illustrate that solar neutrino data do not give 
practical limits. 

For the smaller (and more realistic) values of the e parameters the situation is somewhat more optimistic. In the 
LMA region, neutrino propagation in the Sun remains adiabatic if such effects are present and thus the electron 
neutrino survival probability is controlled by a single parameter, cos29i. The best potential to observe such effects 
with solar neutrinos then seems to be in the E^, ^ few MeV region where 2BiE takes the value cos 26*12, allowing en 
term to be more 'visible'. We illustrate this behavior in Fig. [9] where the percentage change of the electron neutrino 
survival probability. 



-P.e^.e(eil =0) 

is plotted. Future experiments such as the SNO+ experiment [M], currently under construction, should be able to 
probe such solar neutrino energies. 



IV. CONCLUSIONS 



Motivated by the recent low-threshold analysis of the Sudbury Neutrino Observatory and asymmetry measurements 
of the BOREXINO Collaboration of the solar neutrino flux, we reexamined matter effects on the neutrino propagation 
in the Sun. Our analysis shows that solar neutrino data constrains the mixing angle ^13 poorly. Both the BOREXINO 
measurement of the day-night asymmetry of the solar neutrinos and the KamLAND measurement of the reactor 
antineutrinos point out to the LMA region of the neutrino parameter space. In this region solar electron neutrino 
survival probability is not very sensitive to the value of ^13 : to obtain a given value of this probability a wide range 
of values of 6*13 can be used, but permitted variations of 6*12 are limited. 

We next revisited the effects of the physics beyond the Standard Model on the neutrino propagation through the 
Sun. We showed that subdominant Standard Model effects (such as inclusion of a non-zero value of 6*13 or loop 
corrections) can mimic the effects of the physics beyond the Standard Model. For example. Standard Model loop 
corrections are very small, but it is not clear that they are smaller than contributions from new physics. We find that 
, at present, solar neutrino data constrain such effects also rather poorly. 
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Appendix A: Standard Model Corrections to the Two-Flavor Evolution 

Neutrino propagation in matter has been extensively investigated in the context of two-flavor mixing. In this 
appendix we give the reduction of the full Standard Model contribution (i.e., three active flavors and loop corrections) 
to two flavors. We use the following parameterization of the neutrino mixing matrix: 



-23-^13- 
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( Cl3 


5136-**^^ 





(^23 


^23) 





1 





— '5'23 


C23 / 




Ci3 



C12 '5'i2 

-^12 C12 I (A.l) 
1 



where Cij = cos 9ij, Sij = sinOij, and 6cp is the CP-violating phase. Following Ref. [37| we introduce the combina- 
tions 

= C0S6'23^'^- sin 6*23 *r, (A.2) 
'f^ = Sin6l23*^ + COs6'23*r, (A.3) 



which yields the MSW evolution equation 




(A.4) 



where 



El \ (V,^ 

H = T13T12 ( i?2 T\^T\^ + ^IjKm -C23^23Kp I . (A.5) 

S3 / V -C23523^rM CI3KM 



In writing Eq. (A.5 1, by dropping a term proportional to the identity is chosen to be zero. The non-zero potentials 
in this equation are 

V^,,{x) - ^GfN,{x) (A.6) 

and the Standard Model loop correction [55] . 



Tism'^dw \'mw J L mw \2 3 / J 



which is quite small for the solar densities. Note that the initial conditions on and are the same as those on 
^'^ and ^'t-: They are initially all zero. We next perform the transformation 

cos6'i3*e + sin6ii3*^ \ 
-sin6li3*e + cos6li3*^/ 




after which Eq. (A.5 I takes the form 




(A.9) 
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where we dropped a term proportional to the identity. In this equation "H is given by 



iy- A2iCos26ii2 + G 
n = I iA2isin26'i2 + i^ 

i sin 2013 (F^e-i^)e*^^^ 



iA2isin26li2 + 



where we introduced the modified matter potential 



A21 cos 26*12 - G 
iGi3 sin 2023 



23 yT^^ 



and the quantities 



where 



Ai 



G = Vr, 



V = Vr COS^ 



2E 



^13, 



2E ' 



Sin^ 013 COS^ 023 



i sin 2013 (l/^e"i?)e 
-iCi3 sin 2023 V 



—iScF 
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(A.IO) 
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K ^ - sin 013 sin 2023 V^^e 



-iSr 



D = COS^ 023V'rM> 



1 sin^ 013 I cos^ ' 



'23 



sm 



^23 



If one sets V^-p, — 0, 8cp = in Eq. (A.IO) one gets 



n 



'if- A21 COS 2012 iA2i sin 2012 
iA2i sin 2012 -iV + A2i cos 2012 



iKsin20 



13 







i^c sin 2013 


i(A3i + A32) + K 



(A.13) 



Setting the quantity sin 20i3 in the off-diagonal terms of Eq. (A.13 1 yields the result given in Eq. (11.14) [T H [TBI 1461, 147 ) . 



A similar, but more complicated result can be obtained starting from Eq. (A.IO) 
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